1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 300,000 iteration;
- 200,000 burn-in;
- 100 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.00 1.01
## LDevsum 1.00 1.00
## dh0 1.76 3.14
## dh1 2.08 4.06
## dh2 1.22 1.76
## dh3 1.01 1.05
## dl0 1.05 1.11
## dl1 1.44 2.32
## dl2 1.10 1.38
## dl3 1.01 1.02
##
## Multivariate psrf
##
## 1.9
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1218.56835 | 22055.6116 |
| DIC3 | 1159.04687 | 22132.8959 |
| PWAIC | 43.50255 | 239.8298 |
| WAIC | 1186.08453 | 22153.7768 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 300,000 iteration;
- 200,000 burn-in;
- 100 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.04
## LDevsum 1.00 1.01
## dh0 1.01 1.03
## dh1 2.20 6.96
## dh2 1.21 1.75
## dh3 1.59 2.98
## dl0 1.02 1.07
## dl1 1.08 1.33
## dl2 1.01 1.03
## dl3 1.04 1.16
##
## Multivariate psrf
##
## 1.95
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1380.31317 | 22187.433 |
| DIC3 | 1254.95288 | 22235.386 |
| PWAIC | 93.24739 | 295.985 |
| WAIC | 1325.63251 | 22267.482 |